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Ultracold atoms in Raman-dressed optical lattices allow for effective momentum-dependent in- 
teractions among single-species fermions originating from short-range s-wave interactions. These 
dressed-state interactions combined with very flat bands encountered in the recently introduced 
optical flux lattices push the Stoner instability towards weaker repulsive interactions, making it 
accessible with current experiments. As a consequence of the coupling between spin and orbital 
degrees of freedom, the magnetic phase features Ising nematic order. 
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Recently, there has been considerable effort to observe 
the Stoner instability to itinerant ferromagnetism with 
ultracold gases [TH5]. So far, this effort has been fruit- 
less and it has been argued that rapid dimer formation 
at the large coupling strength that is required for fer- 
romagnetism precludes the formation of magnetic do- 
mains [6j[7]. Here we show that atoms subjected to opti- 
cal lattices involving coherent Raman coupling of internal 
states can have a strongly enhanced Stoner instability. 
The ferromagnetic phase appears at much weaker cou- 
pling strength where the gas is less susceptible to dimer 
formation. Furthermore, our results display several in- 
triguing novel phenomena, such as interaction-induced 
phase transitions between distinct Fermi surface topolo- 
gies and nematic ordering, allowing close parallels be- 
tween the physics of cold gases and phenomena in di- 
verse systems such as high temperature superconductors, 
ruthenates and quantum Hall systems [8]. 

Central to our studies are the novel effects that arise 
when atoms are subjected to Raman dressing. Raman 
dressing has recently been used in experiments to create 
artificial gauge potentials [9HTTj and to induce effective 
higher partial-wave interactions among identical bosons 
from short-ranged s-wave interactions [12 . We consider 
an atomic Fermi gas subjected to an optical flux lat- 
tice [13] , in which both of these effects are important. 
The orbital effects of the gauge field cause the lowest en- 
ergy band of the optical flux lattice to be very narrow 
in energy even for a shallow lattice far from the tight- 
binding limit. Interactions among fermionic atoms in this 
lowest band remain sizeable [14]. Simple s-wave interac- 
tions between distinct bare fermions give rise to effec- 
tive interactions of non-zero range among single species 
fermions of the lowest band p~4j [15] . We show that these 
interactions, within this narrow band, cause a ferromag- 
netic transition at a much smaller coupling than in the 
continuum. In view of the coupling of spin and orbital 
motion through the Raman dressing, the ferromagnetic 
transition appears as a change in Fermi surface topol- 
ogy. Furthermore, it is accompanied by a reduction of 
the (spatial) crystal symmetry so also involves nematic 
order [8 . As we describe, this coupling of spin and or- 
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FIG. 1: (Color online) Bandwidth W of the lowest band as a 
function of lattice depth of the optical flux lattice discussed 
in the main text (with — 0.3 and e — 0.4). The inset shows 
the dispersion of the lowest two bands for Vo/Er = 2.25. As 
the lattice is ramped up, the dispersion of the lowest band de- 
velops a minimum at the center of the Brillouin zone causing 
reconstruction of the Fermi surface for non-interacting par- 
ticles. The locations where these reconstructions occur for 
filling v — 1/4 are marked by the dashed lines. For vanishing 
lattice depth Vb = the Fermi circles of the spin-t (red) and 
spin-| (blue) are displaced from each other to the corners of 
the Brillouin zone. 



bital motion allows the magnetic/nematic ordering to be 
readily measured in experiment by band-mapping tech- 
niques. 

We consider the implementation of an optical flux lat- 
tice of Ref. [14], which involves two-photon dressing of 
hyperfine states. We focus on an atomic species with 
groundstate angular momentum F — 1/2, but our key 
ideas are readily extended to atoms with larger F as 
described at the end of the paper. We further restrict 
attention to a quasi-2D geometry, assuming an in-plane 
confinement energy that is large compared to all other 
energy scales. 
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The optical flux lattice is formed by interference of 
three linearly polarized in-plane laser beams with wave 
vectors K\ = — tt/2(-\/3, 1), ^2 = ft/2(V3, — 1) and k<$ = 
ft(0, 1). A fourth circular polarized laser beam, oriented 
perpendicular to the 2D plane, provides the other fre- 
quency required for two-photon Raman coupling of the 
Zeeman-split hyperfine states. 

This optical potential leads to a single particle Hamil- 
tonian 

H = ^i + \4 c (r)i + <x.B(r), (1) 
in which the atom experiences a scalar potential 

V sc = V (3cos 2 (6) - 1) J2 C0 < K J ' r ) ( 2 ) 

3 

and its spin & couples to an effective magnetic field 

B Z = V3V sin 2 (0) sin K' • r ) ( 3 ) 

3 

B x + iB y = eV cos(6)Y,e~ inj ' r ' (4) 

3 

Vq denotes the lattice depth, is the polarization angle of 
the in-plane beams with respect to the surface normal, e 
is proportional to the ratio between Raman coupling and 
scalar potential, and k[ = K\ — K2, k' 2 = ^3 — K\ and 
k' 3 = K2 — K3 [II]- The geometry of the Raman beams 
is such that conversion from spin-^ to spin-^ involves a 
momentum exchange of a^i, or K3. This causes the 
spin character of the Bloch states to vary with crystal 
momentum within the Brillouin zone. (The reciprocal 
lattice basis vectors can be taken to be Gi = K\ — K3 
and G2 = ^2 — ^3.) Notably, for vanishing lattice depth 
Vb = 0, when the energy eigenstates are simply plane 
waves for spin-| and spin-|, the crystal momentum for 
the spin- 4- (t) state with zero kinetic energy is simply Ki 
( — ^2)? or an y equivalent point related by the addition of 
reciprocal lattice vectors. Hence, the crystal momenta of 
the two spin states are displaced from each other within 
the Brillouin zone. In this limit Vb = 0, an unpolarized 
state of non-interacting fermions therefore appears as two 
filled Fermi circles, centered on K\ (—^2) for spin-^ (^) 
and shown in blue (red) in Fig. [I] 

The width of the lowest-energy band is shown in Fig. [I] 
as a function of the overall lattice depth Vb for fixed val- 
ues of e and 6. The bandwidth passes through a min- 
imum at Vo/Er ~ 2, with the recoil energy defined by 
E R = h 2 K 2 /(2m). This is the regime where the optical 
flux lattice best mimics the orbital effects of a uniform 
magnetic field. The lowest energy band is similar to a 
Landau level: with small bandwidth and Chern number 
of one [14]. In the vicinity of this point the positions of 
the band minima change within the Brillouin zone. (Sim- 
ilar features are found in tight-binding models of Chern 



insulators when next nearest-neighbor hoppings are in- 
cluded [15].) This reconstruction is illustrated in Fig. [I] 
by the non-interacting Fermi surfaces shown for a band 
filling of v = 1/4. Note that at Vo = the Fermi surface 
consists of two disconnected circles: these are the spin- up 
and spin-down Fermi surfaces, displaced in crystal mo- 
mentum as described above. 

For non-zero Vo/Er the spin composition of the Bloch 
state continuously varies with crystal momentum. As 
a result, s- wave interactions between spin- up and spin- 
down components lead to effective momentum-dependent 
interactions between fermions in this band. It is remark- 
able that even though we started with a model of short- 
ranged interactions, we obtain an effective theory of in- 
teracting spin- less fermions [13] . 

Atoms restricted to states in the lowest band are de- 
scribed by the effective Hamiltonian 

H lb = 6 k C k C k + ^ ^k 1 k 2 k 3 k4Ck 1 C k 2 C k 3 Ck 4 (5) 

k kik 2 k 3 k4 

where ek is the band dispersion, c k ^ are the fermionic 
field operators for state of crystal momentum k, and 

Vk 1 k 2 k 3 k 4 =#2D / d 2 r ^0 kia (r)0 k2 ^(r)0 k3 ^(r)0 k4a (r) 

J a 

is the effective interaction in the lowest band in terms of 
the eigenfunctions (0kt( r )> ^k|(r)) T of the single particle 
Hamiltonian ([!]). In the following we characterize the 
bare interaction strength by the dimensionless coupling 
parameter g = mgzv/fr 2 . 

To study the effects of interactions, we perform a 
Hartree-Fock (HF) variational approximation which re- 
sults in the energy functional 

E l n k\ = Y 6knk + 2 ^ ^kk^ k n k / (6) 

k kk' 

with Vkk 7 = Vkk'k'k — Vkk'kk 7 - For our zero-temperature 
results, we find the ground states that minimize Eq. (|6| 
for a fixed total number of particles N = ^ k n k . This is 
achieved by setting the occupation numbers n k equal to 
unity for the N orbitals with lowest HF energies 

6c = e k + Y ^kk'^k' • (7) 
k' 

We determine these energies self-consistently by numeri- 
cal iteration, discretizing momenta of the Brillouin zone 
on a fine grid [17 . For our results at non-zero tempera- 
ture T, we instead find n k by minimizing the thermody- 
namic potential Q = E—TS—fiN where \i is the chemical 
potential and S = —ks ^k nk hi(n k ) + (l — n k ) ln(l — n k ) 
is the entropy. The grand canonical potential is station- 
ary (5Q/5nk = 0) when 

nk = e&-M)/^T + I ( 8 ) 



which we numerically iterate to self-consistency with 
Eq. ([7]). We start from different trial states and then 
compare their grand canonical potentials in order to find 
the minimum free energy configuration. 

Our results show a robust ferromagnetic phase for a 
wide range of parameters. To characterize this phase we 
use the magnetization per particle along the z-direction 
as an order parameter 

m z = - [ d 2 r [n t (r)-n ; (r)] . (9) 

In Fig. [2] we show the (modulus of the) calculated or- 
der parameter at zero temperature for a band filling of 
v — 1/4, as a function of 2D coupling g and of lattice 
depth Vo/Er. These results are representative of other 
fillings with v < 1/2. In the absence of any lattice, 
Vo/Er = 0, there is a transition from paramagnet to 
ferromagnet at g c = 2ir. This is the conventional Stoner 
instability for fermions with contact repulsion in 2D [18] . 
This transition appears as a reconstruction of the Fermi 
surface from two Fermi circles to one Fermi circle. As 
the lattice depth is increased, our results show a steady 
decrease of the coupling at which the ferromagnetic tran- 
sition occurs. The minimum coupling for ferromagnetism 
arises for Vq/Er ~ 2, close to the condition for minimum 
bandwidth, Fig. [T] We find that the interaction strength 
g c required is reduced from its free-space value by a factor 
of about 1/4. This is one of the main results of this paper: 
the reduction of coupling as compared to the free-space 
case means that ferromagnetism can be achieved without 
requiring as close an approach to a Feshbach resonance. 
A change of coupling by a factor of 1/4 is expected to 
have a very dramatic reduction in the rate of dimer for- 
mation [EM]. 

In addition to transitions between distinct Fermi sur- 
face topologies, phases of interacting fermions can also 
spontaneously break lattice symmetries. For our model, 
we find that the appearance of ferromagnetism, with non- 
zero magnetization, is also accompanied by a breaking of 
rotational symmetry. We argue that this is a general 
feature of Raman-dressed atomic systems. While spin- 
rotational invariance is explicitly broken by the coupling 
of spin and orbital degrees of freedom, the optical flux 
lattice retains a discrete symmetry. It is invariant under 
a spin- flip combined with a 2tt/6 rotation in real space 

Ue = v x R2*/6- (10) 

The Stoner ferromagnetic transition (effectively a 
Pomeranchuk instability in the spin-channel) causes 
spontaneous symmetry breaking of the Cq spin-rotation 
symmetry arising from Uq down to a residual C3 sym- 
metry associated with Uq = R^/z- This phase transi- 
tion is analogous to the lattice symmetry breaking in the 
electronic Ising nematic phases in solid state materials 
[8, 22-24 . Since the order parameter of the symmetry 
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FIG. 2: (Color online) Magnetization of the ground state of 
interacting fermions in the lowest band of the optical flux 
lattice at uniform filling factor v — 1/4 as a function of lattice 
depth Vo/Er and dimensionless coupling strength g. The 
gray shaded areas on the insets illustrate occupied states in 
the first Brillouin zone. White lines mark transitions between 
unmagnetized states of different Fermi surface topology. 

broken phase is in the 2D Ising universality class, we ex- 
pect this phase to survive to non-zero temperature. 

The phase diagram at non-zero temperature is shown 
in Fig. [3] as a function of chemical potential for a band 
with Vo/Er = 2 and g = 1.9. This shows that the ferro- 
magnetic phase is a robust phase across a range of densi- 
ties and temperatures. The maximum transition temper- 
ature of ksT ~ 0.14W at /i — e m i n ~ 1.54VK corresponds 
to an entropy per particle of S/N ~ 1.1. Entropies of 
this order are being reached in current optical lattice ex- 
periments [25] . 

At the lowest temperatures, we find evidence for the 
appearance of competing phases involving small scale re- 
constructions of the shape (s) of the individual Fermi sur- 
face (s) which can be viewed as Pomeranchuk instabili- 
ties. We do not pursue these here, as we expect these 
low temperature phases to compete with strongly corre- 
lated phases, related to fractional quantum Hall states 
(at specific fillings), which cannot be accessed in HF the- 
ory, but which can commonly co-exist with ferromag- 
netism [28-30 . While our mean- field theory neglects cor- 
relations [26 which can lead to quantitative changes in 
the location of the Stoner instability, detailed numerical 
studies of related models show a robust ferromagnetic 
phase [4j [5j [27] . We expect the reduction of the critical 
coupling strength that we predict to be a robust feature, 
since the shallow optical flux lattice leads to a significant 
reduction in the bandwidth with small decrease in the 
interaction matrix elements. 

Experimental studies of the ferromagnetic transition 




FIG. 3: (Color online) Left: Magnetization for V /E R = 2, 
# = 0.3, e = 0.4 and ^ = 1.9 as a function of temperature and 
chemical potential (scale same for m z as in Fig. [2}. W is the 
bandwidth, and \i — e m in is the chemical potential measured 
from the bottom of the lowest band. The red dots lie on the 
spinodal where the symmetric state becomes unstable (the 
line is a guide to the eye). Right: Occupation numbers rik 
within the first Brillouin-zone for (/x — e m i n )/W = 0.89, 1.69, 
2.19 corresponding to (a), (b), (c). 



we predict will require the use of an atomic species for 
which both strong interactions and Raman coupling can 
be achieved without significant heating. For a hyper- 
fine groundstate of F = 1/2, the natural candidates are 
171 Yb or 199 Hg [14] . For 171 Yb, s-wave contact inter- 
actions between the two states of the lowest hyperfine 
manifold (as spin-^/i) can be conveniently tuned via an 
optical Feshbach resonance [3TJ [32] . Another possibility 
is to use an effective two-level system formed by exploit- 
ing the quadratic Zeeman effect to Raman couple two of 
the hyperfine states in the F = 9/2 groundstate manifold 
of 40 K with neighboring mp [TO] - The interstate inter- 
actions can be conveniently tuned by one of the set of 
magnetic Feshbach resonances that exist for these levels. 

The most direct way to measure the order parameter, 
in Figs. [2] and [3] is by individually imaging the total spin 
populations N^-^. An important practical consequence 
of the coupling of spin and orbital degrees of freedom 
is that total magnetization is not conserved. Hence this 
allows the formation of a macroscopic net magnetization 
starting from an initially unpolarized gas. Measurements 
of the net magnetization will give the average properties 
of the inhomogeneous cloud in the trap. If in addition 
one can measure the local in-situ spin populations n^(r) 
along a contour of fixed filling factor one could then 
map out the phase diagram at different fillings and lat- 
tice depth analogous to recent studies of two-dimensional 
Bose gases [33] . 

Another complementary probe sensitive to the (trap 
averaged) Fermi surface of the dressed lowest band 
fermions is the adiabatic band mapping technique [34] 
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FIG. 4: (Color online) Experimental signatures of adiabatic 
band mapping, (a) Contour plot of the dispersion ek for 
Vo/Er = 1.8, indicating a rhombus shaped reciprocal lat- 
tice unit cell. The Bloch states inside the upper red (lower 
blue) triangle in (a) map to spin-t (40 states with free space 
momentum q illustrated in (b) . For a completely filled lowest 
band one would observe a fully occupied hexagram (b). When 
the atoms pass subsequently through a Stern-Gerlach filter, 
the two spin states can be separately resolved. This would al- 
low clear signatures of the transition from (c) unmagnetized 
to (d) magnetized phase. 



[35] . This probe has the additional advantage that it 
also allows the detection of unmagnetized phases with 
different Fermi surface topologies. Here the lattice po- 
tential is ramped down at a rate slow compared to the 
band gap and fast compared to many-particle dynamics. 
Then the Raman-dressed Bloch states are adiabatically 
mapped onto free-particle plane-wave states of definite 
spin, which can be imaged after time-of- flight expansion. 

The form of this mapping can be deduced by recall- 
ing that for vanishing lattice depth, Vb = 0, the spin-^ 
(t) free-particle state with zero kinetic energy has crystal 
momentum m (—^2) (or any point related by reciprocal 
lattice vectors). Under adiabatic band-mapping, a Bloch 
state whose crystal momentum k is closer to K\ than to 
— K2 is mapped to the spin-^ state with free-space mo- 
mentum q = k — K\\ if k is closer to — k<i the Bloch state 
is mapped to the spin-| state with free-space momentum 
q = k + k,2- This construction is illustrated in Fig. [I] 

For a completely filled lowest band the occupied states 
appear then as a hexagram (superimposed triangles for 
spin ^ and \) after band mapping as shown in Fig. |4jb) . 
A spin resolved image of the cloud of atoms after band 
mapping, obtained either by appropriate detuning of 
imaging lasers or a Stern-Gerlach filter, allows therefore 
for reconstruction of the occupation numbers n^. Finally, 
we note that signatures of the change in Fermi surface 
topology will also appear in the Hall response [36], since 
the Berry curvature of the lowest band of the optical flux 
lattice is non-uniform [37 , so the Hall coefficient is sen- 
sitive to the distribution n^. 
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